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Abstract: We investigate the dynamics of a pair of coincident D5 branes in the 
background of k NS5 branes. It has been proposed by Kutasov that the system with 
a single probing D-brane moving radially in this background is dual to the tachyonic 
DBI action for a non-BPS Dp brane. We extend this proposal to the non-abelian case 
and find that the duality still holds provided one promotes the radial direction to a 
matrix valued field associated with a non-abelian geometric tachyon and a particular 
parametrization for the transverse scalar fields is chosen. The equations of motion 
of a pair of coincident D5 branes moving in the NS5 background are determined. 
Analytic and numerical solutions for the pair are found in certain simplified cases 
in which the U(2) symmetry is broken to U(l) x U(l) corresponding to a small 
transverse separation of the pair. For certain range of parameters these solutions 
describe periodic motion of the centre of mass of the pair 'bouncing off' a finite sized 
throat whose minimum size is limited by the D5 branes separation. 
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1. Introduction 

Kutasovfl, 2] has presented intriguing links between systems of unstable D-branes 
and the DBI effective action of open string tachyon modes of non BPS D-branes. The 
former can be considered for example as a probe BPS D-brane moving in a back- 
ground geometry which breaks all remaining supersymmetry. An example of such a 
geometric background is that due to k coincident NS5 branes [3]. It then emerges 
that one can associate the radial motion of the probe brane in this background with 
that of the open string tachyon of non BPS D-branes. Such an association gives the 
former a geometrical interpretation, hence the notion of 'geometric tachyons'. 

Kutasov's original model was further investigated and extended in [4, 5, 6, 9]. 
Cosmological applications of geometrical tachyons were considered in [10] following 
Sen's original rolling tachyons ideas in [11]. 

In this article we want to investigate what happens if we consider not just a 
single probe D-brane but rather a coincident pair of probe D-branes moving in the 
background of k coincident NS5 branes. For k large, this coincident pair of branes 
can still be regarded as probes in the sense that one may neglect the backreaction 
on the geometry to first approximation. 

This paper is structured as follows: in Section 2 and Section 3 we consider 
different ansatz for the scalar fields which realise the map of the unstable D-brane 
system to one described by a non-abelian tachyonic mode. In section 4 we show the 
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importance of a careful choice of definition for the harmonic function H describing 
the NS5 branes background and we stress the differences between the matrix and 
function approach. In section 5 we make use of a symmetry breaking ansatz to expose 
a simplified version of the non-abelian system and give solutions for the equations 
of motion of the tachyonic field, these will be shown to reduce to the known single 
brane results in the abelian limit. Finally, in Section 6 we provide a small discussion 
of the results and suggest routes for further work. 



2. Multiple D-branes in the NS5-brane background 

Consider a stack of k parallel NS5 branes in type II string theory, stretched in 
the directions x M = (t, x 1 , . . . , x 5 ), // = 0, ... 5, and localised in x m = (x 6 , . . . , x 9 ), 
m = 6, . . . 9. The background fields around k parallel NS5 branes are 

ds 2 = G AB dx A dx B = r]^ v dx^dx u + 5 mn H(x m )dx m dx n (2.1) 
e 2($-$ ) = H(x m ) 

H mnp = -el np d g <S> (2.2) 

where the index A = (/i, m). The function H(x n ) is the harmonic function describing 
k five-branes, H mnp is the field strength of the Kalb-Ramond -B-field and $ is, as 
usual, the dilaton field. For coincident NS5 branes the harmonic function H(x n ) 
reduces to 

kl 2 

H=l + ^ (2.3) 

where r = |x| is the radial coordinate away from the five-branes in the transverse R 4 
labeled by (a; 6 , • • • , x 9 ) and l s = \fa' is the string length. 

We are interested in the dynamics of two coincident BPS D5 branes in the 
background of the five-branes. We can label the world-volume coordinates of the 
D-branes by and by using reparametrization invariance on their world- volume we 
set £ M = x^ . 

The low-energy dynamics of the D5-brane pair is described by a non-abelian 
U(2) gauge theory [12] (see also [13]). The dynamics of the open string sector light- 
est degrees of freedom, namely the adjoint valued scalar fields (X 6 (£ M ), . . . ,X 9 (£ M )) 
which describe the position of the pair in the transverse directions (x 6 , . . . ,x 9 ), the 
non-abelian gauge field as well as the lightest degrees of freedom of the closed 
string sector, namely the metric Gab, the dilaton and the Kalb-Ramond field Bab 
is governed by the non-abelian DBI action 

S = 

- T 5 I d 6 xSTr f e -(*-*°) ^J-det (P + E^ m (Q^ - 5) mn E rw ] + XF^) det(Q%) 
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with 



A = 2vrZ 2 , E AB = G AB + B AB , and Q m n = 6™ + i\[X m , X k }E kn . (2.4) 

The field strength F^ u = d IM A u — d u A^ — i[A^,A v ], P denotes the pullback to the 
brane world-volume and STr denotes the symmetrised trace. 



3. Fuzzy-sphere ansatz for the bulk scalars 

'Fuzzy sphere' configurations for the adjoint scalars X m in the previous non-abelian 
DBI action have been considered in the past [12, 14] . Let us generalise for the 
moment and consider the case of N coincident D5-branes rather than just two and 
consider the following 'fuzzy sphere' ansatz for the transverse scalar fields: 

X* = -R(a;>V , i = 1, 2, 3, (3.1) 

where a 1 give some TV x N matrix representation of the SU(2) algebra 

[a h oij] = 2ie ijk a h . (3.2) 

We define the physical radius of the 5-dimensional transverse space as 

R 2 ( x ») = ^J2 Tr [ X \^Y] = A 2 Ci?(x^) 2 (3.3) 
i=i 

where C is the Casimir of the particular representation of the generators under 
consideration, defined by the identity 



J2^^ i = Ct N (3.4) 



i=i 



Now, given this ansatz the DBI action becomes 



S = -T 5 I cfx STr ( ^=\A + \ 2 Hd a Rd a Ra i a i \Jl + AX 2 R 4 H 2 a i a i ) (3.5) 



with 

hi 2 

H = 1 + ^ (3.6) 
R 2 

where it is understood that H is a function of the physical radius R 2 , and not a 
matrix (this choice is investigated further in Section 4). Note that the symmetrized 
trace in the action ensures that one cannot simply replace all a % a % by the Casimir 
C, there will be ordering issues which spoil this. 
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This action resembles a modified DBI action in flat background of iV non-BPS 
D5-branes proposed in [15], namely 



Sdrt — —T? 



-det (rj ab + Xd a Ti (Q-% d b T^ ) (3.7) 



where 



Qij = l N 8ij - i [T h Tj] (3.8) 

with T\ = To\ and T2 = T02 and there is no sum over i, j. Now we proceed to study 
two limits corresponding to regions of R space where the harmonic function H takes 
two limiting forms, which correspond to the probe D-branes being close and far from 
the NS5's respectively. As will be shown, these limits need to be taken carefully in 
order to preserve forms of the DBI action which have sensible expandable forms. 

3.1 Large radius 

In this section we are looking for a limit of R space in which H — > 1. Hence we need 
R 2 » kl 2 but, to obtain an expandable DBI action which is crucial to performing 
calculations involving the Symmetrised Trace we must also have XR 2 = 2tt1 2 R 2 << 1. 
When this limit is taken the DBI action (3.5) becomes: 

S = -T 5 J (fxSTr (^7= V 7 ! + X 2 d a Rd a Ra rn a m ^l + 4\ 2 R A a n a n ^j . (3.9) 

Expanding both square roots and performing the symmetrised trace manually 
we find that this action is dual to 

S=-T 5 J d 6 xSTr (v(TiTi)Vl + 2T 4 \J I + ~\d a Td a Taia^j (3.10) 

under the map T 4 = 2\ 2 R A CN (where A = and C is the Casimir of the 

N dimensional representation generated by the ai) up to order X 2 R A which given the 
choice of limit means this term is small and higher order terms are progressively less 
important. This is simply the large 1 T expansion of (3.7), where to carry out this 
expansion one needs to be careful in the choice of i = 1, 2 as detailed in [15]. 

The potential takes the form: 

l r (T 2 ) = l--^^ (3.11) 

which is simply the long range gravitational attraction between the D-branes and 
the five-branes. 

Marge here is a slight misnomer, with the map used we still have T 4 << 1, which is important 
for the latter action to be expandable 
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3.2 Small radius 



In this section we are looking for a region of i?-space where 

kl 2 

H ~ (3.12) 

which is achieved for R 2 < kl 2 . However we still want to have a DBI action which 
is expandable, hence we also need R 2 > X 2 kl 2 , which is a sensible enough region 
provided R 2 is not too small compared to kl 2 originally (recalling that A = 27rZ 2 ). 
In this limit the DBI action becomes: 



S = -T 5 f dPxSTr I -^-Wl + A 2 5a a Wa ! «Vl + 4A 2 (fc/ 2 ) 2 aW 
J \vklg\ R 2 



If we set 



R 

T = y/kl s \n—=- (3.13) 
Vkl s 



the previous action becomes, 

S = -T 5 J rf 6 xSTr [e^Ts y/l + 4A 2 (fc/ 2 ) 2 aW \/\ + \ 2 d a Td a T ai c 
which is the tachyon-DBI action with tachyon potential which is corrected from 

T 

the usual by terms which are derived by expanding the action, taking the 

symmetrised trace and matching terms order by order. 



4. Commutative Ansatz 

Here we shall consider a different ansatz to that in the previous section. Inspired 
by [1], where purely radial fluctuations of the fields on the branes give a geometrical 
description of a dual tachyonic system, we re-write the non-abelian action in terms 
of a radial "direction" defined as X m X m = R 2 , and we parametrize the scalar fields 

as 

X m = f m (8, ( j ) ,x)R (4.1) 

where f m are angular functions with f m f m — 1 an d R is an adjoint valued U(2) 
matrix which we rewrite as a linear combination of U(2) adjoint matrices a a in the 
following way 

R = R a (Oa a (4.2) 

where we have also included the U(l) field Rq and defined cto = l^- With this 
parametrization it is clear that the commutator of the scalar fields vanishes [X m , X n ] = 
0, in particular, Q™ = 5™ . Thus in contrast to the fuzzy sphere ansatz of the pre- 
vious section, one might call this a 'commutative' ansatz. 
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The action of the D5-brane pair becomes: 
S = -T 5 J d 6 xSTr ^-±=yj-det + HD^R a D v R h a a a b + F~^)\ (4.3) 

where 

H=l+ -^V (4.4) 
X m X 

where again it is understood that H is an N x N matrix, and 

X m X m = (R a ) 2 + 2R Q R i a a i + (R^f (4.5) 

This action resembles that of two non-BPS D5-branes proposed in [17] in the case of 
vanishing transverse scalar fields: 



S = -T 5 J (fx STr (v(T) yj -deb {r) a p + XD^TD V T + XF^J . 
4.1 Small radii limit 

In the limit in which R ~ Ri ~ R 2 ~ R3 ~ the action (4.3) reduces to 
S = 



(4.6) 



' d 6 xSTr 



. 12 d^R dKR a 2 + 2d il R Q d'"R i a i + d^R^R^m 
1 + — — n _ 73 72 



\/H y (i?o«o) 2 + 2R R l a ai + (RiCki) 

If we set 

T = >/jfeZ a In (iT ot m ) = In R (4.7) 

the above action becomes 

S = - J d 6 x STr (V(T) ^1 + d^Td^ (4.8) 

and in this limit the potential is given by 

V(T) = -^ e ^~s (4.9) 
Vkl s 

If we define T = T m a m then in order to obtain an explicit expression for the different 
components T m of the tachyon matrix we would expand (4.7) and match order by or- 
der each T m components on the l.h.s. with the respective a m component on the r.h.s. 



There is an important point which must be noted here. The map 4.7 is non- 
linear and hence to show that the duality truly holds one needs to show that this 
non-linearity is consistent in the symmetrisation procedure. The result quoted above 
for the dual action is true only under a symmetrisation with respect to the original 
field R, and not the new field T, which is what it would have to be in order for it 
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to be the Tachyon DBI action. In order to show that this action is dual even under 
the symmetrisation procedure one needs to expand the original action in terms of R, 
perform the symmetrisation and then match order by order under a linear map for 
T. For this case this is a hard task due to the difficult powers of aiCHj appearing in 
the expansion and hence we will not show this and aim to investigate it further in 
future work. 



Note that this form of the map allows one to map the fully non-abelian actions, 
including the covariant derivatives. In particular, using eq. (4.7), under the STr we 
have that 

D a R = d a R + i[A a , R] (4.10) 
= Vkl s e^d a T + iVkl s e^[A a ,T] = Vkl s e^D a T (4.11) 

where in the second line we used the usual fact that [f(R),a a ] = f'(R)[R,a a ] for 
f(R) a continuous power series function of a matrix R = R a a a . This means that 

1 \fkl T 

—D a RD a R = D a TD a T + ^[exp D a T}. (4.12) 

The symmetrized trace STr in the action will ensure that the commutator vanishes 
everywhere, so the non-abelian map including the covariant derivatives is realised in 
this limit. 

4.2 Large radii limit 

In the other case, namely when R ~ Ri ~ R 2 ~ R3 — > 00 the map is realized if we 
set 

T = R (4.13) 

and it is trivial to map the components of T with those of R. In this case the tachyon 
potential becomes 

which is the long-range gravitational attraction between multiple D5-branes and the 
NS5 branes. 

In this limit it is trivial to map the covariant derivatives in the action, one 
simply has D a T = D a R and also note that there are now no symmetrisation issues 
in matching the actions. 

4.3 General solution 

Given the ansatz in eq. (4.1), we would like to show that one can find a general 
map for all values of R between the two actions (4.3) and (4.6). In [1], it was shown 
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that for the case of a single probing D-brane (where now R is a function rather than 
a matrix) one could map the two systems by finding an analytical solution to the 
following differential equation: 

% = Vm) (4.i5) 

and by identifying the tachyon potential with the harmonic function H as follows: 

vm ^ ■ (416) 

In the small and large R limits the map gave useful insight into the dynamics of the 
probing brane and provided useful information regarding rolling tachyonic solutions 
[1] and the nature of unstable D-brane systems. In the non-abelian case, the general 
requirement to realise the map is 

STr (HD^RD V R) = STr (\D^TD U T) (4.17) 

When the system is promoted to a non-abelian one such a map is still possible. 
However, one needs to be careful with the choice of definition of H . One possibility 
is that H can be thought of as a matrix, in which case H(RR) depends in a general 
way on the matrix product of R, or we could understand H to depend on R via the 
non-abelian distance H(Tr(R 2 )) so that if is a function and not a matrix. We will see 
through the rest of the paper that the choice is important. Careful string scattering 
calculations should reveal the true form of the Harmonic function appearing in the 
non-abelian DBI action and we think that once these calculations are performed the 
functional form of H will be obtained. However, being unaware of this result in the 
present literature, we decided to pursue both routes and obtain significantly differing 
results. 

4.3.1 H function 

Consider the case where H is chosen to be a function. We will show here that 
analytical solutions for T(R a ) still exist and furthermore that they yield the expected 
single brane results of [1] in the abelian limit. We consider the simplifying case where 
the gauge fields are turned off. 

In this case we define a physical radius as 

3 1 1 

R 2 = J2 ^ xmx ' m = ^ Tr ^ 2 = R 2 + Ri + Ri + R 2 3 (4.18) 

m=l 



With this choice one has 



2 



kll kll 



H <™ = 1 + ^ = 1 + Rl + Rl + Rl + Rl (419) 



2 Notice that another ansatz which makes H a function is H(X m X m ) ~ TrH{X m X m ). This 
ansatz would lead to different results from those we find below and we do not pursue this approach 
any further. 
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In this case we can solve the full map analytically. For every value of R we need to 
solve 

dfj,T = a/ H(R 2 )d fl R . (4.20) 
If we write T = T m a m then for each m = 0, ... 3 we have to solve 



d,T m 



1 + 



M 2 S 



;d„R r > 



R + R 1 + i? 2 + -R3 
In the abelian case with a single D5-brane we would find the solution 



(4.21) 



/ f^§ iR = vwrw + i >„ ffgg^; = i-w (4,2) 

where T kut refers to the tachyon field of the single probe brane case of [1] . By contrast, 
in the non abelian case we have to solve, for example, for the m = component 




Ml 



R 2 + d 2 



dRo 



i^/d 2 + M 2 E isinh" 1 



Rq 

d 



d 2 



d 2 + M 2 



(4.23) 



where we define d 2 = R\ + R\ + R\ for simplicity and E(z, u) is the incomplete elliptic 
integral of the second kind. Although it is trivial to take the limit in which d — > 
on the l.h.s., one has to take care with this limit on the r.h.s. due to divergences 
appearing in the argument of the elliptic integral . In order to explore the differences 
between the abelian and non-abelian case it is instructive to expand the explicit 
expression for the integrand on the l.h.s. in the limit in which R\ + R\ + R\ <ti R 2 , 
Then we obtain: 



T = T kut (R 



+ O 




^sinh" 1 (ag. 



the second term here denotes the non-abelian corrections to the abelian result. 
For reference we write below the full solution for all m 



T 



c(R 



i\l M 2 + R 2 - R 2 m E 



i sinh 1 



R n 



R 2 -Rl 



R 2 -Rl 



M 2 + R 2 -R 2 



for m — 0, 1, 2, 3. where R 2 = R 2 , + R\ + R\ + R\ and R m corresponds to the com- 
ponent of T m one wants to solve for and c(R^ m ) is an integration constant. 
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4.3.2 H Matrix 

In this case we would like to solve the map (4.17) where H is in general a non- 
diagonal matrix. The map is non-trivial (in the case where no particular limit for 
R is taken) unless H is diagonalised, this can be achieved by choosing an a priori 
diagonal ansatz for R. In this case the full Z7" (2) symmetry of the problem would be 
broken to 17(1) <g> U(l). To illustrate this take R = R°a + R 3 a 3 , then 




0\ fdR° + dR 3 



u V (R°-R 3 ) 2 / 

and substituting for R + = R°+R 3 , = R°-R 3 , T+ = T°+T 3 and T_ = T°-T 3 
then one arrives at the map 



kl 2 

dT + = x jl ^dR+ (4.25) 

R+ 



kl 2 

dT_ = x jl-^dR_ (4.26) 

which has as solutions two copies of the solution found in [1]. In particular, the 
action 

S = -T 5 J d 4 xSTr-^= (yi - Hd a Rd a R} (4.27) 

becomes 

S T = -T 5 j d 4 x (e^k y/i - d a T + d a T + + y/l - d a T^d a T^j (4.28) 

which is the U(l) ® U(l) symmetric double copy of the single brane case. This is 
to be expected from a diagonal ansatz, the D-brane probes effectively separate and 
have independent single probe dynamics. 

We believe that general results for arbitrary components of H can be found 
but we will not pursue this any further as the calculations become very involved. 
However, as shown above one can indeed find maps for the H matrix case by taking 
the limits of large/small R first. 



5. Dynamics 

In the case where H is regarded as a diagonal matrix we have seen the effective action 
is just the direct sum of two independent actions each describing the dynamics of 
a single probe D5-brane, which has already been investigated in [1]. Regarding H 
as a function of the non-abelian distance defined in eq. (4.1) produces a dynamical 
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system where there is a non-trivial interaction between the probe branes if we choose 
to separate them (which breaks U(2) — > U(l) x U(l)) . Such an interaction vanishes 
in the flat space limit, as one would expect because then the probe branes are fully 
BPS and no force exists between them whether separated or coincident. 

We take a symmetrical parametrization of the scalar fields, and demand that 
they depend only on time t via 

X m (t) = f m (6,<j ) ,x)R(t) . (5.1) 
Starting from the following action 3 

S = -T 5 J rf 6 xSTr ^=yj-det + Hd^Rd u R)^j , (5.2) 

we make a diagonal ansatz for the scalar field R 

R = R a + R 3 a 3 (5.3) 

and finally we set 

(f) = R0 + R3 

X = Ro-Rs (5.4) 
One finds the action (5.2) reduces to 



S = T 5 I d 6 x-L ( \Jl- #0 2 + ^l-Hx 2 ) (5.5) 



where the harmonic function H is now given by 



H = 1 + 7^ (5-6) 



The equation of motion that follows in the limit in which ~ \ ^ kl 2 is 

<p{xx + <pU-kl 2 J 







(5.7) 

with an analogous one in which \ and are interchanged. Exact solutions to these 
equations are hard to find but one can consider the conservation of the energy which 
results in a simpler first order differential equation. The energy E of the system is 
defined as 

E = P^ + P x x-C (5.8) 



3 We set again the gauge fields to zero manually 
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and we investigate the following ansatz 4 

<P= l -(R Q + C) (5.9) 

X =l(Ro-C) (5.10) 

where C is a constant. In the small Rq limit the conservation of the energy gives 

^ = 2(C + fl8)_47g(C + ^ 
kP, £ 2 («f) 

By imposing reality of the solution one obtains an important inequality 

2kll 42? 

> - 1, (5-12) 



Rl + C 2 ~ E 2 

we see that there are solutions at a critical energy E crit = 2T 5 which can escape to 
infinity. 

The energy equation (5.11) has analytical solutions for C non-zero 



Ro = izCJacobiSN 




rp2 / rp2 97^2^2 

1 + 2 1 -^CHti\-^ + ^C 2 c 1 



E 2 V E 2 —E 2 + 2T 2 C 2 



(5.13) 

where c\ is an integration constant. In Figure 1 we present a plot of this solution 
(with given choice of signs ) for certain values of the parameters E, k, l s , c\ and the 
parameter C = 0,0.01,0.1. They all correspond to the regime where the throat 
approximation to H is valid. The case C = corresponds to the abelian case where 
Ro(t) describes motion which is isomorphic to that of a single probe brane in an 
infinite throat, with energy less than the critical energy required to escape to infinity 
[1]. What is particularly interesting in the case where C ^ is that the solutions 
appear to bounce, at least if we identify the solutions with negative values of Ro as 
separated probe branes moving up the throat. Looking at the harmonic function H 
it is clear that in the case C ^ 0, the geometry seen by the probes is one of a finite 
cutoff throat, with C acting as a cutoff parameter. So the resulting centre of mass 
dynamics of the separated probe pair is equivalent to a single probe brane moving 
in a cutoff throat background. 

In this interpretation, the sub-critical energy probe falls down the throat but 
then reflects off the boundary and back up the throat reaching a certain maximum 
distance, the motion being repeated forever. It's clear from the plots that the period 
of oscillation increases with decreasing C. This makes sense as in the limit C — > we 
recover the solution found in [1] which does not oscillate (at least not in coordinate 

4 by a slight abuse of notation we have used the notation Rq though it is not strictly exactly the 
same as the quantity occurring in eq. (5.3). 
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time t ) but corresponds to an infalling probe brane taking infinite coordinate time 
to reach the throat bottom. By patching the two solutions (which differ by a minus 
sign) together in the regions where Ro(t) is negative one finds an explicit change of 
sign in the velocity Ro(t) of the branes as they reach the throat cutoff, as is expected 
from a perfectly elastic bouncing solution. We now investigate the behaviour of the 
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10 20 30 40 50 

t 

Figure 1: Plots of Rq vs t. The bold curve has C = and corresponds to the abelian 
case. The continuous regular curve corresponds to C = 0.01, and finally the dashed curve 
to C = 0.1. In all cases we have chosen -gf- = 10 , c\ = and kl s = 1 

string coupling with time using the relation 

e 2 * = g s H(Tr(RR)) (5.14) 

In Figure 2 we show a plot of the effective string coupling using the solution (5.13), 
valid in the throat approximation. The thick curve corresponds to the abelian case 
(7 = and shows, as one expects, a rapidly increasing effective string coupling as the 
probe falls down the infinite throat. Thus after some time t = t max the solution is 
no longer within the perturbative string approximation. As argued in [1], the value 
of t max depends analytically on the energy E of the probe and there is an energy 
'window' g s Ts <C E <C T5 for which the probe brane moves in the throat and remains 
within perturbation theory. 

By contrast , the case where C^0 (regular and dashed curves in Figure 2), we 
see the effective coupling as oscillating in time as the probes oscillate in the throat. 
By choosing the value of E and/or C it is possible to control the motion such that 
the string coupling is always in the perturbative regime and for the probes to remain 
in the throat region for all time. 
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Due to the complexity of the solution (5.13) one cannot derive a simple expression 
for a bound on the energy and/or C in order for the above to hold, even for small 
C. Instead one has to use the full expression for the JacobiSN function for C ^ 
and thus we are limited to numerical plots as in Figure 2. 




0.0 ' — 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 L 

10 20 30 40 50 



t 

Figure 2: Plots of the effective string coupling vs t. The bold curve has C = and 
corresponds to the abelian case. The continuous regular curve corresponds to C = 0.01, 
and finally the dashed curve to C = 0.1. In all plots we have chosen -^f- = 10 , c\ = and 
kl 2 s = 1. The value of g s = 0.0001. 

6. Discussion 

In this article we have attempted to generalize the notion of Kutasov's geometric 
interpretation of the open string tachyon, in the scenario where a D-brane is moving 
in a background geometry of k NS5 branes that render the system non-BPS [1]. 
The generalization we investigated considered a pair of coincident probe D5-branes 
moving in this background instead of a single probe D5 brane discussed in [1]. The 
single real geometric tachyon field that appears in the single probe case is, in the 
simplest scenario, related to purely transverse radial motion of the probe. This 
system is abelian in that there is a U{1) gauge theory on the probe brane world 
volume. 

When we consider the case where, for example, one has as a probe two coincident 
D5-branes, then the situation becomes more subtle. Firstly the probe world-volume 
now supports non-abelian U(2) gauge fields and secondly, as is well known, the 
coordinates transverse to this probe stack become matrix valued. This latter phe- 
nomenon raises the question of how one interprets the geometrical quantities such as 
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the harmonic function H sourced by the NS5 branes. In one interpretation, we can 
define a notion of non-abelian distance in the transverse matrix geometry via the 
quantity Tr(X m X m ) where X m are the matrix valued transverse coordinates. Then 
H(X m ) can be thought of as a function via H = H(Tr(X m X m )). Another possible 
interpretation is that H becomes a matrix through its dependence on X m . 

Both definitions seem to give rise to well defined actions since ultimately the 
lagrangians are matrix valued objects in each case and Str is taken over all free 
gauge indices. However as we have shown, the resulting definition of the matrix 
valued geometric tachyon field and the resulting dynamics is different in the two 
interpretations. 

As an illustration of this we saw that in the case where H is treated as a function 
of non-abelian distance defined above, the tachyon map can be found exactly and 
in the limit where the U(2) adjoint valued radial coordinate R is dominated by the 
terms proportional to the 2x2 identity matrix, we recovered the single probe brane 
tachyon map of Kutasov. 

On the other hand, a general solution for the tachyon map in the case of H 
being a matrix is very complicated and its explicit form is not known. However we 
found that at least in the symmetry breaking case where U(2) — > U(l) x U(l) the 
system reduces to two non-interacting copies of single geometrical tachyon fields. By 
contrast, the same U(2) breaking configuration of the probe stack, in the case where 
if is a function and not a matrix, yields a dynamical system involving two coupled 
geometric tachyon fields. 

In this case we found analytic expressions for homogeneous time dependent solu- 
tions at least in the situation where we consider only diagonal degrees of freedom in 
the non-abelian tachyon field, which corresponds to U(2) symmetry breaking. Inter- 
estingly we found oscillating or 'bouncing' solutions in this case where the separation 
parameter between the two D5 probes acting as an effective cutoff on the NS5 brane 
infinite throat. 

It would be very interesting to find (even numerically) dynamical solutions which 
involve the full non-abelian degrees of freedom in the U(2) valued tachyon field in 
the action eq (5.2) including non-vanishing gauge fields. Another extension could be 
to look at different arrangements of background NS5 branes other than the point like 
ones considered in this paper. For example one can also consider the k NS5 branes 
arranged around a ring of finite radius. This is a known supergravity solution and 
the corresponding metric and harmonic function are known [18] . This would extend 
to the non abelian case the results found in [5]. 
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